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The sensitivity achievable by a pair of VIRGO detectors to stochastic and isotropic gravitational 
wave backgrounds of cosmological origin is discussed in view of the development of a second VIRGO 
interferometer. We describe a semi-analytical technique allowing to compute the signal-to-noise 
ratio for (monotonic or non-monotonic) logarithmic energy spectra of relic gravitons of arbitrary 
slope. We apply our results to the case of two correlated and coaligned VIRGO detectors and 
we compute their achievable sensitivities. The maximization of the overlap reduction function is 
discussed. We focus our attention on a class of models whose expected sensitivity is more promising, 
namely the case of string cosmological gravitons. We perform our calculations both for the case of 
minimal string cosmological scenario and in the case of a non-minimal scenario where a long dilaton 
0^ ■ dominated phase is present prior to the onset of the ordinary radiation dominated phase. In this 

ON ' framework, we study possible improvements of the achievable sensitivities by selective reduction of 

ON ■ the thermal contributions (pendulum and pendulum's internal modes) to the noise power spectra 

of the detectors. Since a reduction of the shot noise does not increase significantly the expected 
sensitivity of a VIRGO pair (in spite of the relative spatial location of the two detectors) our findings 
support the experimental efforts directed towards a substantial reduction of thermal noise. 
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I. THE PROBLEM AND ITS MOTIVATIONS 



It is well known that every variation of the background geometry produces graviton pairs which are stochastically 
distributed and whose logarithmic energy spectra represent a faithful snapshot of the (time) evolution of the curvature 
scale at very early times |jj . Indeed, one of the peculiar features of stochastic graviton backgrounds is that their energy 
spectra extend over a huge interval of (present) frequencies. This feature can be appreciated by comparing the graviton 
backgrounds with other backgrounds of electromagnetic origin (like the cosmic microwave background [CMB] ) . The 
Q\ I analysis of the CMB background (together with its spatial anisotropies) is relevant for very large (length) scales Q 
ON ■ (roughly ranging between the present horizon [i.e. 10 -18 Hz] and the horizon at decoupling [i.e. 10~ 16 Hz]). Since 
' gravitational interactions are much weaker than electromagnetic interactions they also decouple much earlier and, 
^j*. therefore, the logarithmic energy spectra of relic gravitons produced by the pumping action of the gravitational field 
can very well extend for (approximately) twenty five orders of magnitude in frequency Q . From the physical point of 
view, this observation implies that the energy spectra of relic gravitons can be extremely relevant in order to probe the 
past history of the Universe in a regime which will never be directly accessible with observations of electromagnetic 
backgrounds. 

In spite of the fact that the nature of the production mechanism is shared by different types of models Jj]] , the specific 
amplitudes of the energy spectra can very well change depending upon the behavior of the background evolution. An 
example in this direction are logarithmic energy spectra increasing in frequency Q]. Different theoretical signals 
(with different spectral distributions) lead to detector outputs of different amplitudes. We are facing a non-linear 
problem where a change in the detector signal can be determined either by an improvement in the features of the 
detector or by a different functional form of the logarithmic energy spectrum Therefore, in order to evaluate the 
performances of a given detector one has to choose the specific functional form of the logarithmic energy spectrum. 
A possible choice is represented by scale invariant spectra Hf7|. Another rather interesting choice is represented b 
tilted ( "blue" S ) spectra whose energetical content is typically concentrated at frequencies larger than the mHz 
String cosmological models |l(J are yet another interesting theoretical laboratory leading usually to sizable theoretical 
signals in the operating window of wide band interferometers (WBI) |nj . A possible detection of these backgrounds 
would represent an interesting test for cosmological models inspired by the low energy string effective action. 

Every measurement in cosmology turns out to be difficult for different and independent reasons. The CMB 
anisotropy experiments have to cope with the mandatory subtraction of different electromagnetic foregrounds which 
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can be much larger than the "cosmological" signal one ought to detect. In order to detect gravitational waves of 
cosmological origin with terrestrial measurements we are facing similar problems. 

The signal induced in the detector output by stochastic gravitational waves backgrounds is indistinguishable from 
the intrinsic noise of the detector itself. This implies that, unless the amplitude of the signal is very large, the only 
chance of direct detection of these backgrounds lies in the analysis of the correlated fluctuations of the outputs of, 
at least, two detectors affected by independent noises. The problem of the optimal processing of the detector outputs 



required for the detection of the stochastic background has been considered by various authors 12 13|] and it was also 
reviewed in ref. p4| . 

Suppose, indeed, that the signal registered at each detector can be written as (we limit ourselves to the case of two 
detectors (i = 1,2)) 

Si = hi(t) + ni(t) , (1.1) 

where we have indicated with n the intrinsic noise of the detector, and with h the gravitational strain due to the 
stochastic background. By assuming that the detector noises are stationary and uncorrelated, the ensemble average 
of their Fourier components satisfies 

«(/)%(/')) = \ &{f-f')5ijS®<\f\) > (1-2) 

where S^fl/I) is usually known as the one-sided noise power spectrum and is expressed in seconds. Starting to the 
signals S\ and s%, a correlation "signal" for an observation time T can be defined in the following way: 

S = [ dt [ dt' Sl {t)s 2 (t')Q(t-t') (1.3) 

J -T/2 J-T/2 

where Q is a filter function that depends only by t — t' because we assume that n and h are both stationary. The 
optimal choice of Q corresponds to the maximization of the signal-to-noise ratio associated to the "signal" S. Under 
the further assumptions that detector noises are Gaussian, much larger in amplitude than the gravitational strain 



and statistically independent on the strain itself, it can be shown 12-14 that the signal-to-noise ratio in a frequency 
range (/ m , / M ) is given by []: 

SNR2 _ 3tf 2 j— I [*" 7 2 (/)^w(/) I' 7 ' ri4 s 

- ^ FVT \L df 7^mr)S • (L4) 

where Hq is the present value of the Hubble parameter and F is a numerical factor depending upon the geometry of 
the two detectors. In the case of the correlation between two interferometers F — 2/5, howev er, in the correlation of 
detectors of different geometry, F ^ 2/5 (see Appendix A for details about this point). In Eq. Jl.4| ), the performances 

(1 2) 

achie vabl e by the pair of detectors are certainly controlled by the noise power spectra (NPS) Sn ■ However in 
Eq. fll.4p , on top of NPS, there are two important quantities. The first one is the theoretical background signal 
defined through the logarithmic energy spectrum (normalized to the critical density p c ) and expressed at the present 
(conformal) timoP] rj 

^Gw(/i1o) = — - ~fr~7 = ^(Va)^(f,Vo) ■ (1.5) 
Pc din/ 

The second one is the overlap reduction function "f(f) |l3| , |l4| ] which is a dimensionless function describing the reduction 
in the sensitivity of the two detectors (at a given frequency /) arising from the fact that the two detectors are not 
in the same place and, in general, not coaligned (for the same location and orientation "/(f) — 1). Since the overlap 



reduction function cuts-off the integrand of Eq. (1.4) at a frequency which approximately corresponds to the inverse 
separation between the two detectors, it may represent a dangerous (but controllable) element in the reduction of the 
sensitivity of a given pair of detectors. 



1 In order to avoid possible confusions we stress that the definition of the SNR is the one discussed in [^| and it is essentially 
the square root of the one discussed in p"^-|l4||. 

2 In most of our equations we drop the dependence of spectral quantities upon the present time since all the quantities 
introduced in this paper are evaluated today. 
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Various ground-based interferometric detectors are presently under construction (GEO LIGO-LA, LIGO-WA 
|l6f , TAMA , VIRGO ) . Among them, the pair consisting of most homogeneous (from the point of view of 
the noise performances) detectors with minimum separation is given by the two LIGOs (VIRGO and GEO are even 
closer, but they have different performances for what concerns the NPS). However, this separation (~ 3000 km) is 
still too large. The overlap reduction function 7(/) for the pair LIGO-LA— LIGO-WA encounters its first zero at 64 
Hz, falling off (swiftly) at higher frequencies, i.e., right in the region where the two LIGOs, at least in their initial 
version, have better noise performances. 

Recently, within the european gravitational wave community, the possibility of building in Europe an interferometric 
detector of dimensions comparable to VIRGO has received close attention [^9|. Therefore, there is a chance that in 
the near future the VIRGO detector, now under construction at Cascina (Pisa) in Italy will be complemented by 
another interferometer of even better performances very close (at a distance d < 1000 km) to it. In this paper 
we examine in detail the possible improvements in the VIRGO sensitivity as a result of direct correlation of two 
VIRGO-like detectors. Furthermore, since technological improvements in the construction of the interferometers can 
be reasonably expected in the next years, it is easy to predict that also VIRGO, as for the LIGO detectors, will 
gradually evolve toward an advanced configuration. For this reason we also examine the possible consequences of a 
selective improvements of the noise characteristics of the two detectors on the obtained results. 

In order to evaluate precisely the performances of a pair of VIRGO detectors we will use the following logic. First 
of all we will pick up a given class of theoretical models which look particularly promising in view of their spectral 
properties in the operating window of the WBI. Secondly we will analyze the signal-to-noise ratios for different regions 
of the parameter space of the model. Finally we will implement some selective reduction of the noises and we will 
compare the results with the ones obtained in the cases where the noises are not reduced. We will repeat the same 
procedure for different classes of models. 

The results and the investigations we are reporting can be applied to spectra of arbitrary functional form. The 
only two requirements we assume will be the continuity of the logarithmic energy spectra (as a function of the present 
frequency) and of their first derivative. We will also give some other examples in this direction. 

In order to make our analysis concrete we will pay particular attention to the evaluation of the performances of a 
pair of VIRGO detectors in the case of string cosmological models | p^ , |il| . 

The plan of our paper is then the following. In Section II we introduce the basic semi-analytical tecnique which 
allows the evaluation of the SNR for a pair of WBI. In Section III we will evaluate the performances of a pair of 
VIRGO detectors in the case of string cosmological models. In Section IV we will show how to implement a selective 
noise reduction and we will investigate the impact of such a reduction in the case of the parameter space of the models 
previously analyzed. Section V contains our final discussion and the basic summary of our results. We collect in the 
Appendices some technical results useful for our analysis and other interesting complements to our investigation. 



II. SNR EVALUATION 



In the operating window of the VI RGO detectors the theoretical signal will be defined through the logarithmic 
energy spectrum reported in Eq. (1.5). In the present Section we shall not make any specific assumption concerning 
u>(f) and our results have general applicability. We will only assume that it is a continuous function of the frequency 
and we will also assume that its first derivative is well defined in the operating window of WBI. This means that u>(f) 
can be, in principle, a non-monotonic function. 



A. Basic Formalism 



The noise power spectrum of the VIRGO detector is well approximated by the analytical fit of Ref . ^ZE , namely 

f<fb 



£»(/) 



So 



3C 



Si 



1 + 



f>fb 



(2.1) 



where 



Ei = 3.46 x 10~ 6 

So = 10~ 44 s , f a = 500Hz , f b = 2Hz , E 2 = 6.60 x 10~ 2 

= 3.24 x 10~ 2 . 
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In order to compute reliably (and beyond naive power counting arguments) the SNR we have to specify the overlap 
reduction function "f(f). The relative location and orientation of the two detectors determines the functional form of 
7(/) which has to be gauged in such a way that the overlap between the two detectors is maximized (i.e. 7(/) — 1 
for most of the operating window of the two VIRGO). Moreover, the two interferometers of the pair should also be 
sufficiently far apart in order to decorrelate the local seismic and electromagnetic noises. Since the precise location 
of the second VIRGO detector has not been specified so far Jl9|, we find useful to elaborate about this point by 
computing the overlap reduction functions corresponding to two coaligned VIRGO interferometers with different 
spatial separations. The results of these calculations are reported in Fig. [|. Needless to say that these choices are 
purely theoretical and are only meant to illustrate the effects of the distance on the performances of the VIRGO pairn. 



FIG. 1. We report the overlap reduction function(s) for the correlation of the VIRGO detector presently under construction 
in Cascina (43.6 N, 10.5 E) with a coaligned interferometer whose (corner) station is located at: A) (43.2 N, 10.9 E), d — 58 
km (Italy); B) (43.6 N, 4.5 E), d = 482.7 km (France); C) (52.3 N, 9.8 E), d = 958.2 km (Germany). The third site (C) 
corresponds to the present location of the GEO detector. Notice that from A to C the position of the first zero of j(f) gets 
shifted in the infra-red. See also Appendix A concerning this last point. 

The curves labeled with A, B, and C shown in Fig. |l| correspond to different distances d between the site of the 
VIRGO detector (presently under construction in Cascina, near Pisa) and the central corner station of a second 
coaligned VIRGO interferometer. Let us now look at the position of the first frequency /j for which 7(/j) = for 
each of the curves. We can notice that by increasing d (i.e., going from A to C) the value of /) gets progressively 
shifted towards lower and lower frequencies, linearly with d. This means that, for the specific purpose of the detection 
of a stochastic gravitational waves background, the position of the first zero of the overlap reduction function cannot 
be ignored. On a general ground we would like /, to be slightly larger than the frequency region where the sensitivity 
of the pair of wide band detectors is maximal. In the explicit examples presented in this paper we will focus our 
attention on the case A. The other two configurations will be the subject of a related investigation |20|. 



3 For illustrative purposes, we assumed that a distance of about 50 km is sufficient to decorrelate local seismic and e.m. noises. 
Such a hypothesis is fair at the present stage and it is certainly justified within the spirit of our exercise. However, at the 
moment, we do not have any indication either against or in favor of our choice. 
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B. SNR versus phenomenological bounds on the graviton spectrum 



By inserting the parametrization (1.5) into Eq. fll.4| ) we can write 



SNR 2 = J^°- VT J > ( 2 - 2 ) 

where we introduced the (dimension-less) integral 

Here the integration variable is v = ///o, with /q a generic frequency scale within the operating window of the 
interferometer, and the integration domain is restricted to the region f m < f < fu (i-c, v m < v < vy{). In the 
following we will choose /o = 100 Hz and, taking into account the frequency behavior of 7(/) (see Fig. Q), we can 



assume /m = 10 kHz (i.e., vyy = 100). The lower extreme / m is put equal to the frequency /& entering Eq. (2.1) 
(i.e., i/ m = 0.02). 
For the chosen values of /o and So (see Eq. (|2.lD) one has: 



^ 4 -° X T 1(r7 r^V /2 SNR 2 . (2.4) 



j v T 

Since we will often refer to this formula we want to stress its physical meaning. Suppose that the functi onal form of 



oj{f) is given. Then the numerical value of the integral J can be precisely computed and, through Eq. (2.4), fl can 
be estimated. This quantity, inserted in Eq. fll.5| ), determines for each frequency / the minimum f^GW detectable 
(for an observation time T, with a signal-to- noise ratio SNR) by the correlation of the two detector outputs. 

In the next section, Q will be compared with two other quantities: Q th and f2 max . The first is the theoretical 
value of the normalization of the spectrum, while the second represents the largest normalization compatible with the 
phenomenological bounds applicable to the stochastic GW backgrounds. These quantities are of different nature and 
in order to be more precise let us consider an example. 

Suppose, for simplicity, that we are dealing with a logarithmic energy spectrum which is a monotonic function of 
the present frequency. Suppose, moreover, that the spectrum decreases sufficiently fast in the infra-red in order to 
be compatible both with the pulsar timing bound and with the CMB anisotropics bounds. Then the most relevant 
bound will come, effectively, from Big-Bang nucleosynthesis (BBN) |2^-^]. Therefore, in this particular case, we will 
have that is determined by demanding that 

hi J G w(/,%)dln/ < 0.2/1^(770) =s 5 x 10~ 6 , (2.5) 

where f2 7 (?7o) = 2.6 x 10~ 5 h^ 2 is the present fraction of critical energy density stored in radiation. According to 
our definition, f2 max is the maximal normalization of the spectrum compatible with the previous inequality, namely, 



^H max * , 1 = din/. (2.6) 

1 JfnB 



Notice that f ns ~ 10~ 10 Hz is the present value of the frequency corresponding to the horizon at the nucleosynthesis 
time; /max stands for the maximal frequency of the spectrum and it depends, in general, upon the specific theoretical 
model. If the spectrum has different slopes, will be determined not only by the nucleosynthesis bound but 

also by the combined action of the CMB anisotropy bound [plpjfl and of the pulsar timing bound ]2q] . Indeed, we 
know that the very small fractional timing error in the arrival times of the millisecond pulsar's pulses implies that 
JIqw ^ 10 -8 for a frequency which is roughly comparable with the inverse of the observation time along which pulsars 
have been monitored (i.e., uj p ~ l/r o b s = 10 -8 Hz). Moreover, the observations of the large scale anisotropies in the 



microwave sky 25 imply that the graviton contribution to the integrated Sachs- Wolfe effect has to be smaller than 
(or at most of the order of) the detected amount of anisotropy. This observation implies that Ogw < 6.9 x 10 -11 
for frequencies ranging between the typical frequency of the present horizon and a frequency thirty of forty times 
larger. In the case of a logarithmic energy density with decreasing slope the Q will be mainly determined by the 
Sachs- Wolfe bound and it will be the maximal normalization of the spectrum compatible with such a bound. 
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On a general ground, we will have that fl * < f2 max , namely the theoretical normalization of the spectrum is 
bounded, from above, by the maximal normalization compatible with all the phenomenological bounds. Therefore, 
the mismatching between these quantities can be interpreted as an effective measure of the theoretical error in the 
determination of the absolute normalization of the spectrum. 

Since u>(f) enters (in a highly non-linear way) into the form of J (as defined in Eq. (|2.3|)), the corresponding 



in Eq. (2.4) will be different for any (specific) frequency dependence in Lu(f). The consequence of this statement is 
that it is not possible to give a general (and simple) relationship between the sensitivity at a given frequency, the 
spectral slope and the (generic) theoretical amplitude of the spectrum. However, given the form of the theoretical 
spectrum, the phenomenological bounds (depending upon the theoretical slope) will fix uniquely the theoretical error 
and the maximal achievable sensitivity. So, if we want to evaluate the performances of the VIRGO pair we should 
pick up a given class of theoretical models (characterized by a specific functional form of u)(f)) and compute the 
corresponding sensitivity. The same procedure should then be repeated for other classes of models and, only at the 
end, the respective sensitivities can be compared. 

III. PRIMORDIAL GRAVITON BACKGROUND VERSUS VIRGO*VIRGO 

We can consider, in principle, logarithmic energy spectra with hypothetical analytical forms and arbitrary normal- 
izations. If the logarithmic energy spectrum is either a flat or a decreasing function of the present frequency JHL we 
can expect, on general grounds, that the theoretical signal will be of the order of (but smaller than) 10~ 15 [£7| for 
present frequencies comparable with the operating window of the VIRGO pair. This happens because of the combined 
action of the Sachs- Wolfe bound together with the spectral behavior of ther infra-red branch of the spectrum produced 
thanks to the matter-radiation transition. Of course this observation holds for models where the graviton production 
occurs because of the adiabatic variation of the background geometry 

In order to have large signals falling in the operating window of the VIRGO pair we should have deviations from scale 
invariance for frequencies larger than few mHz. Moreover, these deviations should go in the direction of increasing 
logarithmic energy spectra. This is what happens in the case of quintessential inflationary models |J. In this case, 
however, as we discussed in a previous analysis |>|, the BBN bound put strong constraints on the theoretical signal 
in the operating window of the VIRGO pair. 

Another class of model leading to a large theoretical signal for frequencies between few Hz and 10 kHz is represented 
by string cosmological models [[i 10 llj]. Therefore, in order to evaluate the performances of the VIRGO pair and in 



order to implement a procedure of selective noise reduction we will use string cosmological spectra. 

A. Minimal models of pre-big-bang 

In string cosmology and, more specifically, in the pre-big-bang scenario, the curvature scale and the dilaton coupling 
are both growing in cosmic time. Therefore the graviton spectra will be increasing in frequency instead of decreasing 
as it happens in ordinary inflationary models. 

In the context of string cosmological scenarios the Universe starts its evolution in a very weakly coupled regime with 
vanishing curvature and dilaton coupling. After a phase of sudden growth of the curvature and of the coupling the 
corrections to the tree level action become important and the Universe enters a true stringy phase which is followed by 
the ordinary radiation dominated phase. It should be stressed that the duration of the stringy phase is not precisely 
known and it could happen that all the physical scales contained within our present Hubble radius crossed the horizon 
during the stringy phase as pointed out in Q |. 

The maximal amplified frequency of the graviton spectrum is given by ]4]JTl[| 

h{ m ) ~ 64.8 VffT f — J GHz (3.1) 

where n r is the effective number of spin degrees of freedom in thermal equilibrium at the end of the stringy phase, 
and gi = M s /Mpy where M s and Mpi are the string and Planck masses, respectively. Notice that g\ is the value 



4 An exception to this assessment is represented by cosmic strings models leading to a flat logarithmic energy spectrum for 
frequencies between 10 -12 Hz and 10~ 8 Hz [ ps| . Another possible exception is given by the gravitational power radiated by 
magnetic (and hypermagnetic) E9(] knot configurations at the electroweak scale ]30|. 
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of the dilaton coupling at the end of the stringy phase, and is typically of the order of 10 -2 -j- 10 _1 |32(| . In order 
to red-shift the maximal amplified frequency of the spectrum from the time r\\ (which marks the beginning of the 
radiation dominated evolution) up to the present time we assumed that the cosmological evolution prior to rjo and 
after 771 is adiabatic. Minimal models of pre-big-bang are the ones where a dilaton dominated phase is followed by a 
stringy phase which terminates at th e on set of the radiation dominated evolution. In the context of minimal models, 
the function ui(f) introduced in Eq. (|l.5| ) can be written as 



"(/) 



-2/3 



J_ 

Is 

3-/3 



where, 



1 + *„ 
fi 

P = 



2/3-3 _ I ln L 



2 fs 



-1 2 



/</. = - 

Zs 



fs < f < fl 



(3.2) 



ln (gi/gg 
In z. 



(3.3) 



In this formula z a = f\ / f s and g s are, respectively, the red-shift during the string phase and the value of the coupling 
constant at the end of the dilaton dominated phase. The first of the two branches appearing in Eq. (3.2) is originated 
by modes leaving the horizon during the dilaton dominated phase and re-entering during the radiation dominated 
phase. The second branch is mainly originated by modes leaving the horizon during the stringy phase and re-entering 
always in the radiation dominated phase. The theoretical normalization 



fi 



th 



2.6 ff 2 



10 3 

n r 



1/3 



^7(^0) , 



(3.4) 



multiplied by u)(f) (as given in Eq. (^2)) leads to the theoretical form of the spectrum. Notice that n r is of the order 
of 10 2 -7- 10 3 (depending upon the specific string model) and it represents a theoretical uncertainty. 

However, as anticipated in the previous section, the theoretical normalization of the spectrum should be contrasted 
with the one saturating the BBN bound (i.e., fi ). This quantity is obtained by Eq. (2.6), where in the case under 
consideration 



T = Id 



ith 



Id = 



U J 



Is = 



fs 



(3.5) 



The analytical expressions of 2d and X s are reported in Appendix B. We have to bear in mind that in the intermediate 
frequency region of the graviton spectra an important bound comes from the pulsar timing measurements. Therefore, 
if one ought to consider rather long stringy phases (i.e., large z s ), the BBN constraint should be supplemented by the 
requirement that fiQw(10~ 8 Hz) < 10~ 8 pq |. We will come back to this point later. 

Following the explicit expression of the function oj(f), Eq. (2.4) can be re-written as follows: 



4 x 10" 



1 yr 



1/2 



SNR 



(3.6) 



where, introduced the following notation 

J k = I dv- 



(2), 



ln fc v 



k = 0,1,2,3,4 



J± 



m(3-2/3) 



Av 



7 2 (/o^) 



.±m(3-20) 



1,2 



(3.7) 



C d = l + zf- 3 + - lni/ a , 

one has 
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z s~ 2/3 ( A 3 2 1 1 A 1/2 

Jd = ' S ^3 yCdJn - 2C d J 1 + -jC d J 2 - -jpdJz + J^J4 

/ \ 1/2 

"l I "l "l "l "l / 



The previous expressions are general in the sense that they are applicable for a generic value of f s . If f m < f s < Jm 
then both J s and Jd give contribution to the sensitivity. If, on the other hand f s < f m (i.e., a long stringy phase) 
the main contribution to the sensitivity will come from J s . The integrals appearing in Jd iS have to be evaluated 
numerically. In all our calcu lations we will assume that both VIRGO detectors are characterized by the same (rescaled) 



NPS (reported in Eq. (gj)) 



The main steps of our calculation are the following. We firstly fix g\ and for each pair (z s , gi/g s ) (within the range 
of their physical value) we compute CI (for T = 1 yr and SNR = 1), and CI . We then compare these two quantities 

-th 



to the theoretical normalization given in Eq. (3.4). If CI will be larger than CI (but smaller than CI ) we will say 
that the theoretical signal will be "visible" by the VIRGO pair. In this way we will identify in the plane (z s , gi/g s ) 
a visibility region according to the sensitivity of the VIRGO pair. The theoretical error on the border of this region 

7T max T\ *h 

can be estimated by substituting \ l to \l . 

To illustrate this point we consider a specific case. The value of the coupling at the end of the stringy phase can be 
estimated to lie between 0.3 and 0.03 p2|. The knowledge of g\ will not fix uniquely the theoretical spectrum which 
does also depend on the number of relativistic degrees of freedom at the end of the stringy phase. Therefore, the 
theoretical error in the determination of the absolute normalization of the spectrum could be also viewed as the error 
affecting the determination of n r . In all the plots shown we will take, when not otherwise stated, g\ = 1/20 and 
n r = 10 3 as fiducial values. Different choices of g\ will lead to similar results. We will also assume that the overlap 
reduction function associated with the pair is the one reported in the curve A of Fig. Jl| 

In Fig. U (top left) we report the result of our calculation for the ratio between CI and CI as a function of gi/g s 
and \ogz s . The contour plot (bottom left) shows the region of the plane (logz s , gi/g s ) where this ratio is greater 
than 1, i.e. the maximal visibility region allowed by the BBN bound. In the opposite case, i.e., CI /CI < 1, the 
VIRGO pair is sensitive to a region excluded by the BBN. In the right part of Fig. || we go one step further and 

we plot the ratio between CI and CI. The shaded area in the contour plot (bottom right) is the region of the plane 

(logz s , gi/g s ) where the conditions Q th /Cl > 1 and Q max /Q > 1 are simultaneously met. The shaded area in this 
plot defines the visibility region of the VIRGO pair assuming the theoretical normalization of the spectrum. From 
Fig. H by ideally subtracting the shaded area of the left contour plot from the shaded area of the right contour plot 
we obtain an estimate of the theoretical error. The results we just presented can be obviously recovered for different 
values of g\ close to one. However, if gi gets too small (and typically below 1/25) the visibility area gets smaller and 
smaller eventually disappearing. 

The visibility regions appearing in Fig. ^ extend from intermediate values of z s (of the order of 10 s ) towards large 
values of z s (of the order of 10 18 ). Notice that for our choice of gi, f s can become as small as 10~ 8 for z s of the order 
of 10 18 . As we recalled in the previous Section, this frequency corresponds to the inverse of the observation time 
along which pulsar signals have been monitored and, therefore, for this frequency, a further "local" bound applies to 
the logarithmic energy spectra of relic gravitons. This bound implies that fiGw(10 _8 Hz) < 10~ 8 . In our examples, 
the compatibility with the BBN bound implies also that the pulsar timing constraint is satisfied. Given our choice 
for gi we can clearly see that the visibility regions depicted in Fig. |2| extend for values of g s which can be as small as 
1/160 (or as small as 1/60 in the case of right part of Fig. ||). 



B. Non-minimal models of pre-big-bang 

In the context of minimal models of pre-big-bang, the end of the stringy phase coincides with the onset of the 
radiation dominated evolution. At the moment of the transition to the radiation dominated phase the dilaton seats 
at its constant value. This means that g\ ~ 0.03 0.3 at the beginning of the radiation dominated evolution. As 
pointed out in |3l| , it is not be impossible to imagine a scenario where the coupling constant is still growing while the 
curvature scale starts decreasing in time. In this type of scenario the stringy phase is followed by a phase where the 
dilaton still increases, or, in other words, the coupling constant is rather small at the moment where the curvature 
starts decreasing so that g% <C 1. 
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og z. 
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FIG. 2. We report the ratios f2 max /f2 (left) and ^2 th /S7 (right) as a function of gi/g 3 and logz s (fi is calculated for T = 1 
yr and SNR = 1). The lower contour plots show the regions where these ratios are greater than 1. The shaded area (bottom 
right) represents the region where the combination of the theoretical parameters is such that the corresponding H th does not 
violate the BBN bound. As we can see the visibility region is reduced. The difference between the shaded area in the right plot 
and the one in the left plot measures the error made by assuming as normalization of the spectrum not the theoretical one but 
the maximal one compatible with the BBN. The value z s — 10 s roughly corresponds to f s ~ fo- Notice that log denotes not 
the Neperian logarithm but the logarithm in ten basis. 

After a transient period (whose precise duration will be fixed by the value of gi), we will have that the radiation 
dominated evolution will take place when the value of the coupling constant will be of order one (i.e., g r ~ 1). 



An interesting feature of this speculation is that the graviton spectra will not necessarily be monotonic |31| (as the 
ones considered in the previous analysis). We then find interesting to apply our considerations also to this case. 
The function u>(f) in the non- minimal model described above is given by pl| 



"(/) 



, 2/V5 / , \ 4 
9r\ (J 

.91 



1 - In 



fs 



2/V3 



2-er 



2+cr 



i -ml 



fr < f < fs = 
fs < f < fl 



Zs 



(3.9) 



where, in the present case 

fi - 64.8 VffT 



l/2v/3 



n r 



1/12 



GHz 



fr 



-2/V3 



f 



1 • 



(3.10) 



The frequency f r corresponds to the onset of the radiation dominated evolution. If we adopt a purely phenomcnological 
approach we can say that f r has to be bounded (from below) since we want the Universe to be radiation dominated 



5 Notice that the form of reported in |3lj differs from our expression only by logarithmic correction whose presence is, 
indeed, not relevant. We kept them only for sake of completeness. 
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not later than the BBN epoch. Thence, we have that f r > f ns . Recalling the value of the nucleosynthesis frequency 
and assuming that g r ~ 1 this condition implies g\ ^,8.2 x 10~ 16 . This simply means that in order not to conflict 
with the correct abundances of the light elements we have to require that the coupling constant should not be too 
small when the curvature starts decreasing. Notice that for frequencies / < f r the spectrum evolves as f~ 3 . The 
ultra-violet branch of the spectrum is mainly originated by modes leaving the horizon during the stringy phase and 
re-entering when the dilaton coupling is still increasing. 

Concerning the non- minimal spectra few comments are in order. Owing to the fact that g± can be as small as 
10 -15 we have that the highest frequency of the spectrum can become substantially smaller than in the minimal 
case. Moreover, the spectrum might also be non-monotonic with a peak at f s . Looking at the analytical form of the 
spectrum we see that this behavior occurs if a > 2. A non-monotonic logarithmic energy spectrum (with a maximum 
falling in the sensitivity region of the VIRGO pair) represents an interesting possibility. 

The results of our calculation for g\ = 10~ 12 , n r = 10 3 , g r = l,and a > 2 are reported in Fig. ||[ As done in the 
case of minimal spectra we analyse the visibility window in the plane of the relevant parameters of the model. As we 
can see from the left part of Fig. ||| the region compatible with the BBN is rather large but it shrinks when we impose 
the theoretical normalization ( right part of Fig. |3j) which is always smaller than the maximal normalization allowed 
by BBN. 




log z, log z, 

FIG. 3. In order to make clear the comparison with the visibility region of the minimal models, we report f2 max /fi (left) and 
Q. /fi (right) as a function of a and of the logz s in the non-minimal scenario. Notice that we took gi = 10 -12 , n r — 10 3 , 
and g r = 1. As for Fig. the shaded areas in the lower contour plots represent the region where each ratio is greater than 1, 
and, in the case of the right plot, also the BBN is satisfied. 

It is interesting to compare directly the three dimensional plots appearing in Fig. || with the corresponding three 
dimensional plots of Fig. 0. We can see that the regions of parameter space where Cl™^ jCl and CI /O are larger than 
one is larger in the case of minimal models. However, the shaded region in the case of minimal models corresponds 

to ratios fi max /fi and CI /CI which can be 3 or 2, respectively. On the other hand the shaded region in the case of 

Fig. H corresponds to ratios Sl max /0 and CI /CI which can be, respectively, as large as 50 or 25. So, in the latter case 
the signal is larger for a smaller region of the parameter space. 

As we stressed in the previous Section, CI represents the sensitivity of the VIRGO pair to a given spectrum whose 
functional form is given by ui(f). One might be interested, in principle, in the sensitivity of the VIRGO pair at a 
specific frequency /*. This can be easily computed by multiplying CI by u(f*). In Fig. ^ we show the sensitivity of 
the VIRGO pair at the frequency /* = 100 Hz, both, for the minimal and non-minimal models considered in the 
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present Section. One can easily discuss the same quantity for any other frequency in the operating window of the 
VIRGO detectors. 




FIG. 4. We report the logarithm of the sensitivity of the VIRGO pair at 100 Hz for T = 1 yr and SNR = 1 in the case of 
minimal (left plot) and non-minimal (right plot) energy spectra. 

IV. NOISE REDUCTION AND THE VISIBILITY REGION OF A VIRGO PAIR 

There are two ways of looking at the calculations reported in this paper. One can look at these ideas from a 
purely theoretical perspective. In this respect we presented a study of the sensitivity of a pair of VIRGO detectors 
to string cosmological gravitons. There is also a second way of looking at our exercise. Let us take at face value 
the results we obtained and let us ask in what way we can enlarge the visibility region of the VIRGO pair. In this 
type of approach the specific form of graviton spectrum is not strictly essential. We could use, in principle, any 
motivated theoretical spectrum. As we stressed, we will use string cosmological spectra because, on one hand, they 
are theoretically motivated and, on the other hand, they give us a signal which could be, in principle detected. Of 
course, there are other well motivated spectra (like the ones provided by ordinary inflationary models). However, the 
signal would be, to begin with, quite small. 

In this Section we will then consider the following problem. Given a pair of VIRGO detectors, we suppose to be 
able, by some means, to reduce, in a selective fashion, the contribution of a specific noise source to the detectors 
output. The question we ought to address is how the visibility regions will be modified with respect to the case 
in which the selective noise reduction is not present. We will study the problem for th e pa ir of VIRGO detectors 
considered in the previous Sections, i.e., for identical detectors with NPS given in Eq. ( |2.l| ), and characterized by 
the overlap reduction function of the case A of Fig. |l|. As for the theoretical graviton spectrum we will focus our 
attention on the case of minimal models considered in Section III. A, with the same parameters used to produce Fig. 
g. Also here, the quantity fl will be computed for T = 1 yr and SNR = 1. 

As shown in Section II the NPS is characterized by three dimension-less numbers £1,2,3, and two frequencies /„ 
and /t,. Roughly, Si and £2 control, respectively, the strength of the pendulum and pendulum's internal modes noise, 
whereas £3 is related to the shot noise (see Ref. |33| for an accurate description of the phenomena responsible of the se 
noises). Below the frequency fb the NPS is dominated by the seismic noise (assumed to be infinitum in Eq. ( |2.l| )). 
The frequency f a is, roughly, where the NPS gets its minimum. The frequency behavior of this three contributions 
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and of the total NPS is shown in Fig. [5|. The stochastic processes associated with each source of noise are assumed 
to be Gaussian and stationary. 
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FIG. 5. The analytical fit of the rescaled noise power spectrum E n defined in Eq. (|2.l[) in the case of the VIRGO detector. 
With the full (thick) line we denote the total NPS. We also report the separated contribution of the three main (Gaussian and 
stationary) sources of noise. 

In the following, without entering in details concerning the actual experimental strategy adopted for the noise 
reduction, we will suppose to be able to reduce each of the coefficients Si by keeping the other fixed. In order to 
make our notation simpler we define a "reduction vector" 



P = (Pl,P2,P3) , 



(3.1) 



whose components define the re duct ion, respectively, of the seismic, thermal and shot noises with respect to their 
fiducial values appearing in Eq. (2.1) (corresponding to the case p = (1,1,1)). 

As shown in Fig. [5] the pendulum noise dominates the sensitivity of the detectors in the low frequency region, 
namely below about 40 Hz. In Fig. || we report the results of our calculation for the case p = (0.1, 1, 1). Here the 
parameters of the theoretical spectrum are exactly the same as in Fig. ^. The only change is given by a reduction 
of the pendulum noise. From the comparison between Fig. ^ and Fig. |^, we see that the visibility region in the 
parameter space of our model gets immediately larger especially towards the region of small g s . This enlargement is 
quite interesting especially in terms of Q, /f2. 

In the frequency region between 50 and 500 Hz the performances of the detectors are, essentially, limited by the 
pendulum's internal modes noise. The results obtained for a selective reduction of this component are summarized 
in Fig. |?j, where the pendulum and shot noises are left unchanged but the internal modes component is reduced by 
a factor of ten (in Fig. |J). As we can see the visibility region gets larger and the increase in the area is comparable 
with the one obtained by selecting only the pendulum noise. 

Finally, for sake of completeness, we want to discuss the case of the shot noise, i.e., the noise characteristic of the 
detector above 500 Hz. Our results for p^ = 0.1 are reported in Fig. ||. As we can see by comparing Figs. ^, [?], and 
8 we gain much more in visibility by reducing the thermal noise components than by reducing the shot noise. In Fig. 
3 the shot noise is reduced by one tenth but the visibility region does not increase by much (left plot). This result is 
consequence of the fact that, as shown by Fig. |5|, the shot noise contribution to the NPS starts to be relevant for / ~ 1 
kHz, i.e., in a frequency region where the overlap between the detectors begins to deteriorate (see Fig. 0). In Figs. 
^ and the thermal noise is reduced by one tenth and the increase in the visibility region is, comparatively, larger. 
This shows, amusingly enough, that a reduction in the shot noise will lead to an effect whose practical relevance is 
already questionable at the level of our analysis. Notice that a selective noise reduction can be also discussed in the 
case of a purely flat spectrum (20| . 
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FIG. 6. We report the ratios f2 max /f2 (left plots), and 1 /fi (right plots) in the case in which the shot noise and the noise 
related to the pendulum's internal modes are not reduced, whereas the pendulum noise is diminished by a factor of ten with 
respect to the values quoted in Eq. (2.1), i.e., p — (0.1, 1, 1). 



a °™ / n 
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FIG. 7. We report the result of selective reduction in the case where the noise cause by the pendulum's internal modes is 
reduced by a factor of ten, whereas the pendulum and shot contributions are left unchanged, i.e., p = (1,0.1, 1). 
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FIG. 8. We report the same quantities discussed in Fig. ^|for the case p = (1, 1,0.1). Thhe shaded areas in the lower plots 
are the relevant visibility regions which should be compared with the shaded regions in the lower plots of Figs. |^ and [?]. By 
direct comparison we can argue that a reduction in the shot noise (by a factor of ten) is not as efficient as a reduction, by the 
same amount, in the thermal noise components. 

In order to conclude this Section we want to show the combined action of the simultaneous reduction of both the 
components of the thermal noise. In Fig. ^, owing to the results of our analysis we kept the shot noise fixed but 
we reduced both the thermal and seismic noises by a factor of ten. Clearly we observe a consistent increase in the 
visibility region. However, even if a combined reduction of these components cannot be achieved we want to stress 
that already a reduction of the pendulum's internal modes noise alone (by one tenth) can be of relevant practical 
interest. 



V. DISCUSSION AND EXECUTIVE SUMMARY 



There are no compelling reasons why one should not consider the appealing theoretical possibility of a second 
VIRGO detector coaligned with the first one. Moreover, recent experimental suggestions seem coherently directed 
towards this goal [[[9). While the location of the second detector is still under debate we presented a theoretical 
analysis of some of the scientific opportunities suggested by this proposal. 

We focused our attention on possible cosmological sources of relic gravitons and we limited our attention to the 
case of stochastic and isotropic background produced by the adiabatic variation of the backgound geometry. In 
the framework of these models we can certainly argue that in order to have a large signal in the frequency window 
covered by VIRGO we have to focus our attention on models where the logarithmic energy spectrum increases at large 
frequencies. Alternatively we have to look for models where the logarithmic energy spectrum exhibits some bump in 
the vicinity of the VIRGO operating window. If the logarithmic energy spectra are decreasing as a function of the 
present frequency (as it happens in ordinary inflationary models) the large scale (CMB) constraints forbid a large 
signal at high frequencies. In the case of string cosmological models the situation seems more rosy and, therefore, we 
use these models as a theoretical laboratory in order to investigate, in a specific model the possible improvements of 
a possible VIRGO pair. The choice of a specific model is, in some sense, mandatory. In fact, owing to the form of 
the SNR we can immediately see that different models lead to different SNR not only because the amplitude of the 
signal differs in different models. Indeed, one can convince himself that two models with the same amplitude at 100 
Hz but different spectral behaviors between 2 Hz and 10 kHz lead to different SNR. 
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FIG. 9. We illustrate the case of a simultaneous reduction of Ei and £2 by a factor 10, whereas E3 is the same of Eq. ( |2,l[ ), 
i.e., p = (0.1,0.1,1). 



In order to analyze the sensitivity of the VIRGO pair we described a semi-analytical technique whose main advantage 
is to produce the sensitivity of the VIRGO pair to a theoretical spectrum of arbitrary slopes and amplitudes. The 
theoretical error is estimated, in our approach, by requiring the compatibility with all the phenomenological bounds 
applicable to the graviton spectra. As an intersting example, we asked what is the sensitivity of a VIRGO pair 
to string cosmological spectra assuming that a second VIRGO detector (coaligned with the first one) is built in a 
european site. By assuming that the second VIRGO detector has the same features of the first one we computed the 
SNR and the related sensitivity achievable after one year of observation in the case of string cosmological spectra. 

By using the string cosmological spectra as a theoretical laboratory we then studied some possible noise reduction. 
Our main goal, in this respect, has been to spot what kind of stationary and stochastic noise should be reduced 
in order to increase the visibility region of the VIRGO pair in the parameter space of the theoretical models under 
considerations. Our main result is that a selective reduction of each of the three main sources of noise is not equivalent. 
A reduction in the shot noise by a factor of ten does not increase significantly the visibility region of the VIRGO pair. A 
selective reduction of the thermal noise components is far more efficient. In particular, we could see that a reduction (of 
one tenth) of the pendulum's internal modes increases the visibility region of four times. The simultaneous reduction 
of the two components of the thermal noise leads to an even more relevant increase. 

The construction of a second VIRGO detector coaligned with the first one and an overall reduction of the thermal 
noise of each detector of the pair leads to what we called "upgraded VIRGO" program. The results presented in this 
paper are obtained in the case of a particularly promising class of theoretical models but can be generally applied to 
any logarithmic energy spectrum with similar qualitative results. However, owing to the non-linearities present in the 
evaluation of the SNR it would not be correct assess that they hold, quantitatively, without change. We hope that 
our results and our suggestions may turn out to be useful in the actual process of design of the upgraded VIRGO 
program p9| . 
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APPENDIX A: THE NORMALIZATION OF THE OVERLAP REDUCTION FUNCTION 



In this Appendix we discuss the reduction in sensitivity due to the fact that, in general, these detectors will 
not be either coin cide nt or coaligned. This effect is quantified by the (dimensionless) overlap reduction function j(f) 
appearing in Eq. ( |l.4| ) . Suppose that we have a gravitational wave propagating along a generic direction characterized, 
in spherical coordinates, by the unit vector Cl = (cos cp sin 9, sin <p sin 9, cos 9) . If we now introduce a pair of orthogonal 
unit vectors directed in the plane perpendicular to Q 

m(£l) = (cos0 cos 9, simp cos6>, — sin#) h((l) = (sin</>, — coscp, 0), (B-l) 

the polarization tensors can be written, in terms of the polarization angle ip of the GW, as 

£ + {Cl,ip) = e + (n)cos2^> - e x (A)sin2V> 

e x (n,tp) = e + (f2)sin2?/> + e x (Q) cos 2tp , (B.2) 

where 

e + (n)=rn(jl)®m(n)-n(0)®n(n) e x (Cl) = m{Cl) ® h{Q) + h{fl) ® m(fi) (B.3) 
with the normalization 

Tr{e A (Cl)e A '{Q)} = 25 AA ' . 
If the graviton background is isotropic and unpolarized we will have that 

A 

where Ar = r\ — f% is the separation vector between the two detector sites, F A is the pattern function characterizing 
the response of the i-th detector (i — 1, 2) to the A = +, x polarization, and the following notation 

f dO <hp . . m K . 

has been introduced to indicate the average over the propagation direction (0, 0) and the polarization angle ip. The 
normalization factor i* 1 is given by: 

F = J2 <F A (f 1 A4>)F A {f 2 A^)>n^2, (B.6) 

A 

where the notation 1 = 2 is a compact way to indicate that the detectors are coincident and coaligned and, if at least 
one of them is an interferometer, the angle between its arms is equal to 7r/2 (L-shaped geometry). In this situation, 
by definition, -f(f) = 1. When the detectors are shifted apart (so there is a phase shift between the signals in the two 
detectors) , or rotated out of coalignment (so the detectors have different sensitivity to the same polarization) it turns 
out that: |7(/)| < 1. 

The pattern functions (or orientation factors) of a GW detector can be written in the following form 

F A (r,n,iP) = Tr{D(f)e A (n,i;)} (B.7) 

where the symmetric, trace-less tensor D(r) describes the orientation and geometry of the detector located at r. 

The tensor D(f) depends upon the geometrical features of the detector. For instance, in the case of an interferometer, 
indicating with u and v the unit vectors in the directions of its arms, one has: 

D(r) = - {u{f) ® u(f) - v{f) ® v{f)} . (B.8) 

In the case of the lowest longitudinal mode of a cylindrical GW antenna with axis in the direction determined by the 
unit vector Z, one has 

D{f) = i{f) (g) 1(f) - \l , (B.9) 



1G 



where I is the unit matrix. Finally, in the case of the lowest five degenerate quadrupole modes (m 
a spherical detector, the corresponding tensors are 



-2) of 



D {+1) (f) = - 1 -g x (f) 
L>(+ 2 ) (f) = ie+(f) D 



1 



2^/3 



{2/+(f) - e+(f)} 



(B.10) 



(-2) m 



■ r ><(f) 



where 



f + (f) — 777(f) <E> rh{f) — f <8> ?' 
<7 + (r) = fi(f) ® h(f) — r ® f 



f x (f) = 777(f) ® f + f ® m(f) 
g x (f) = n(f ) ® f + r (g) ft(f) , 

and e +,x (f) are the tensors of Eq. ( |B.3j ) written in terms of the unit vectors 777(f) and ft(f) lying on the plane 
perpendicular to f . From these expressions for the tensors D lJ and interpreting each of the five modes of a sphere as 
a single detector, it is possible to show that in the case of coincident detectors one has: 



<Fi 4 (f,fi»^ 2 s (f,f2,^) >^= c 12 S 



AB 



(A,B 



(B.ll) 



where C12 depends only on the geometry and the relative orientations of the two detectors. The corresponding values of 
F (see Eq. (B.6)) for the three different geometries considered (interferometer, cylindrical bar, sphere) are summarized 
in Table |[ By introducing the following notation 



Ar = ds 



2wfd, 



where s is the unit vector along the direction connecting the two detectors and d is the distance between them, it can 
be shown [Ol that the overlap reduction function assumes the following form (Dk = D(fk)): 



7(/) = Po{6) D{ 3 D 2ij + pi{5) D? D 2 % Sj s k + p 2 (r?) D{ 3 £>*' s t s 3 s k Sl 



where 



Pa 
Pi 
Pi 



(5) = -U 
v ' F5 2 



2S 2 
-46 2 
S 2 



-AS 



2 




jo 




20 




h 




35 




h 





(B.12) 



(B.13) 



with jk{S) the standard spherical Bessel functions: 



Jo(S) 



sin 6 



Ji(S) 



jo(5) - cos 5 



.n{S) = 3— j Q (6). 



TABLE I. The normalization factor F for three different geometries of the detectors: interferometer (ITF), cylindrical bar 
(BAR), and sphere (SPH). A * denotes entries that can be obtained from the symmetry of the table. 
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APPENDIX B: BBN BOUNDS 



In the case of minimal models the integrals determining the analytical expression of the BBN bound are given by: 



+ ^ + (f ) 3 [2tf + 6* + 18) 

-6z s (z s +6) m ^ +9z 2 ln 2 ^ 

Js Js 



o 2/3-6 -2/3 

3 2/3(3-2/?) 2/3-6 2/3 1 ' ' 

In the case of non-minimal models the integrals determining the BBN bound are given by 

Ji = A(a, z s ) + B(a,z s ) lnz s + C(a, z s ) In 2 z s , 

1i = ^ (<- 2 + ^ ) - V 4 ) (1 + ln^) 2 , (B.2) 

where and z T — f\/f r and 

A ^ z ') = ~ 16 (,72 _ 4)3 { 13 ^ - 4) 3 - 4 z; 4 (a + 2) 3 (2a 2 - 10a + 13) 

+ 4z; 4(1+a) (a - 2) 3 (2cr 2 + 10a + 13) - z- 2a (13a 6 - 172a 4 + 832a 2 - 1664) I , 

B ^ *•) = 4^-4)2 { 2 + 2 ) 2 ^ - 5) - 2 z; Aa (a 2) 2 (2a + 5) - 5 zf (a 2 - 4) 2 } , (B.3) 

4-2ct 

Zs) = 2 (a 2 -4) { 2 " ^ ( a - 2) + a ^ (a 2 - 4))} . 
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